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We review constraints on additional Z' bosons predicted in supcrsymmetric (SUSY) E$ 
models from electroweak experiments — Z-pole experiments, m w measurements and the 
low-energy neutral current (LENC) experiments. Four representative models - \, ip, rj, v 
models — are studied in some detail. We find that the improved data of parity violation 
in cesium atom, which is 2.2-cr away from the Standard Model (SM) prediction, could 
be explained by the exchange of the heavy mass eigenstate Z2 in the intermediate state. 
The improvement over the SM can be found in Xj *7: v models, where the total x 2 of 
the fit to the 26 data points decreases by about five units, owing to the better fit to 
the atomic parity violation. Impacts of the kinetic mixing between the U(l)y and 
U(l)' gauge bosons on the x 2_ana ly s i s are studied. We find that the Z 1 model with 
(/3e, <5) = (— 7r/4, 0.2), where 0e is the mixing angle between Z x and Z^ bosons and 
8 denotes the kinetic mixing, shows the most excellent fit to the data: the total x 2 
decreases by about seven units as compared to the SM. We introduce the effective mixing 
parameter f, a combination of the mass and the kinetic mixing parameters. The 95% 
CL lower mass bound of Z2 can be shown as a function of £. A theoretical prediction on 
( and the U(l)' gauge coupling gE is studied for the Xi^PyV an d v models by assuming 
the minimal particle content of the SUSY Eq models. 

1. Introduction 

The presence of an additional Z' boson is predicted in a certain class of grand 
unified theories (GUT) with a gauge group whose rank is higher than that of the 
Standard Model (SM). The supersymmetric (SUSY) Eg models are the promising 
candidates which predict the additional Z'-boson at the weak scale 0. Because 
E 6 is a rank-six group, it can have two extra U(l) factors besides the SM gauge 
group. A superposition of the two extra U(l) groups may survive as the U(l)' gauge 
symmetry at the GUT scale. The U(l)' symmetry may break spontaneously at the 
weak scale through the radiative corrections to the mass term of the SM singlet 
scalar field □. 

In general, the additional U(l)' gauge boson Z' can mix with the hypercharge 
U(l)y gauge boson through the kinetic term at above the electroweak scale, and also 
it can mix with the SM Z boson after the electroweak symmetry is spontaneously 
broken. Through those mixings, the Z' boson can affect the electroweak observables 
at the iT-pole and the lU-boson mass m w . Both the Z-Z' mixing and the direct Z' 
contribution can affect neutral current experiments off the Z pole. The presence of 



2 Looking for Z' bosons in Supper symmetric E§ Models through Electroweak Precision Data 



an additional Z' boson can be explored directly at pp collider experiments. 

In this review article, we report constraints on Z' bosons in the SUSY Eq mod- 
els from electroweak experiments based on the formalism in Refs. ^[ ^|. Constraints 
on tha Z^. bosons from electroweak experiments have been studied by several au- 
thors Q'130'ffil. Especially a special attention has been paid to this subject t£l after 
the new analysis of parity violation in cesium atom has led to the improved data 
of the weak charge QwQ^Cs) Ej> which is 2.2-a away from the SM prediction (see 
Table ^). The analysis given in this article updates their studies by allowing for 
an arbitrary kinetic mixing liini3'L3 between the Z' boson and the hypercharge B 
boson. The constraints on the Z' bosons can be found by using the results of Z-pole 
experiments at LEP1 and SLC, and the m w measurements at Tevatron and LEP2. 
Also the low-energy neutral current (LENC) experiments lepton-quark, lepton- 
lepton scattering experiments and atomic parity violation (APV) measurements - 
constrain the direct exchange of Z' boson. 

It has been found ill that the lower mass limit of the heavier mass eigenstate Z-i 
is obtained as a function of the effective Z-Z' mixing term £, which is a combination 
of the mass and kinetic mixings. In principle, £ is calculable, together with the gauge 
coupling <?£, once the particle spectrum of the E§ model is specified. We show the 
theoretical prediction for Q and in the SUSY Eg models by assuming the minimal 
particle content which satisfies the anomaly free condition and the gauge coupling 
unification 3 . 

This paper is organized as follows. In the next section, we review the additional 
Z' boson in the SUSY Eg models and the generic feature of Z-Z' mixing in order to 
fix our notation. We show that the effects of Z-Z' mixing and direct Z' boson con- 
tribution are parametrized by the following three terms: (i) a tree- level contribution 
to the T parameter 113, T new , (ii) the effective Z-Z' mass mixing angle £ and (iii) a 
contact term ff^/c^m^ which appears in the low-energy processes. In Sec. 3, we 
collect the data of electroweak experiments which will be used in our analysis. We 
also present the theoretical framework to calculate the electroweak observables. In 
Sec. 4, we show constraints on the Z' bosons from the electroweak data. The pres- 
ence of non-zero kinetic mixing between the U(l)y and U(l)' gauge bosons modifies 
the couplings between the Z' boson and the SM fermions. We discuss impacts of the 
kinetic mixing term on the x 2 -analysis. The 95% CL lower mass limit of the heavier 
mass eigenstate Z 2 in four representative models - x> tp, 77, v models - is given as a 
function of the effective Z-Z' mixing parameter £. The ^-independent constraints 
from the low-energy experiments and those from the direct search experiments at 
Tevatron are also discussed. In Sec. 5, we find the theoretical prediction for £ in 
X, ip, T), v models by assuming the minimal particle contents. Stringent Zi boson 
mass bounds are found for most models. Sec. 6 is devoted to summarize this paper. 

2. Z-Z' mixing in supersymmetric E§ model 



a Consequence in the case of the maximal particle con tent which preserve the perturbative unifi- 
cation of the gauge couplings has been studied in Ref. Il4l 
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2.1. Z' boson in super symmetric E$ model 



Since the rank of Eq is six, it has two U(l) factors besides the SM gauge group 
which arise from the following decompositions: 



E 6 D SO(10) x U(l)^ 

DSU(5)xU(l) x xU(l) f 



(1) 



An additional Z' boson in the electroweak scale can be parametrized as a linear 
combination of the U(l)^ gauge boson Z^ and the U(i) x gauge boson Z x as \3 



Z' = Z x cos Pe + Zi/> sin (3e- 
In this paper, the following Z' models are studied in some detail: 



(3e 





tt/2 




tan~ 1 (\/l5) 


model 


X 


i> 


V 


V 



(2) 



(3) 



In the SUSY-i?6 models, each generation of the SM quarks and leptons is embedded 
into a 27 representation. In Table 1, we show all the matter fields contained in a 
27 and their classification in SO(10) and SU(5). The U(l)' charge assignment on 
the matter fields for each model is also given in the same table. The normalization 
of the U(l)' charge follows that of the hypercharge. Besides the SM quarks and 
leptons, there are two SM singlets v c and S, a pair of weak doublets H u and Hd, 
a pair of color triplets D and D in each generation. The rj model arises when Eq 
breaks into a rank-5 group directly in a specific compactification of the heterotic 
string theory 0. In the v model, the right-handed neutrinos v c are gauge singlet 
and can have large Majorana masses to realize the see-saw mechanism E3. 

The U(l)' symmetry breaking occurs if the scalar component of the SM singlet 
field develops the vacuum expectation value (VEV). It can be achieved at near 
the weak scale via radiative corrections to the mass term of the SM singlet scalar 
field. Recent studies of the radiative U(l)' symmetry breaking can be found, e.g., 
in Ref. |. 

Several problems may arise in the Eq models from view of low-energy phe- 
nomenology. For example, the presence of the baryon number violating operators 
give rise to too fast proton decay, or the absence of the Majorana neutrino mass 
terms (except for the v model) requires a fine-tuning of the Dirac neutrino mass in 
order to satisfy experimentally observed neutrino mass relations. Some approaches 
to these problems are summarized in Ref. [lj. In the following, we assume that these 
problems are solved by an unknown mechanism. Moreover we assume that all the 
super-partners of the SM particles and the exotic matters do not affect the radiative 
corrections to the electroweak observables significantly i.e., they are assumed to be 
heavy enough to decouple from the weak boson mass scale. 



2.2. Phenomenological consequences of Z-Z' mixing 
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Table 1. The hypercharge Y and the U(l)' charge Qe of all the matter fields in a 27 for the 
X,ip,V and v models. The classification of the fields in the SO(10) and the SU(5) groups is also 
shown. The value of U(l)' charge follows the hypercharge normalization. 

SO(10) SU(5) field Y 2\/6Q x ^/WJbQ^ ~Q U 




1 v c -5 +1 -| 



10 5 H u +i +2 -2 +| -yf\ 




If the SM Higgs field carries a non-trivial U(l)' charge, its VEV induces the Z-Z' 
mass mixing. On the other hand, the kinetic mixing between the hypercharge gauge 
boson B and the U(l)' gauge boson Z' can occur through the quantum effects below 
the GUT scale. After the electroweak symmetry is broken, the effective Lagrangiap 
for the neutral gauge bosons in the SU(2) L x U(l)y x U(l)' theory is given by tHl 

gauge — ^ uv 2~ ~ 4 ^ 

+m% z ,Z»Z^ + l -m z Z»Z u + \m\,Z'»Z^ (4) 

where F^ v [F = Z, Z' , A° , B) represents the gauge field strength. The Z-Z' mass 
mixing and the kinetic mixing are characterized by rn 2 zz , and sin^, respectively. In 
this basis, the interaction Lagrangian for the neutral current process is given as 

Cnc = -Y.{z QfJ^-ffaK + azToTi* 1 {1% - Qf a sin 2 9w) f a Z u 

f,a 

+g E Q f £T a rf a z' l ) , (5) 
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where gz = g/ cos9w — 9Y / sitlOw- The U(l)' gauge coupling constant is denoted 
by gE in the hypercharge normalization. The symbol f a denotes the quarks or 
lcptons with the chirality a (a — L or R). The third component of the weak 
isospin, the electric charge and the U(l)' charge of f a are given by ij , Qf a and 
, respectively. The U(l)' charge of the quarks and leptons listed in Table 1 
should be read as 



QT 



Q S e (f = e,u,d). (6) 



~ ^E 

The mass eigenstates (Z\, Zi, A) is obtained by the following transformation; 



Z \ / cos £ + sin £ sin 9w tan % 
Z' J = sin£/cosx 
A J \ — sin £ cos 8w tan \ 



sin £ + cos £ sin Ow tan \ 
cos £ / cos x 
— cos £ cos Qw tan x 



Zi 

Z-2 

A 

(7) 
(8) 

The physical 



= m z( c i + S£Swt x ) +m z A — +2m zz , — (c$ + s^swt x ), (9a) 



Here the mixing angle £ is given by 
tan2£ = - 



2c x (m zz , 



+ s w s x m 2 z ) 



• n z> - ( c l - s w s V) m z + 2s w s x m zz> 
with the short-hand notation, c x = cosx, s x — sinx and sw = sin#vF 
masses m Zl and m Z2 (m Zl < m Z2 ) are given as follows; 



m| 2 = m|(cgs w i x - sj) 2 + m z , I — ) + 2m zz ,^- (c e s w t x - sj), (9b) 

\ c x / c x 

where = cos£, = sin£ and t x — tan%. The lighter mass eigenstate Z\ should 
be identified with the observed Z boson at LEP1 or SLC. The excellent agreement 
between the current experimental results and the SM predictions at the quantum 
level implies that the mixing angle £ has to be small. In the limit of small £, the 
interaction Lagrangians for the processes Z\ % fafa are expressed as 



{l 3 fL ~Q fa sm 2 9, 



If a — Qf a sin 2 



faZ\fi 

9e . 



f a Z; 



where the effective mixing angle £ in Eq. ( 10a ) is given as 

9E 



(10a) 
(10b) 

(11) 



gz cosx 

In Eq. (pp|), the effective U(l)' charge Q^* is introduced as a combination of 
and the hypercharge Yf a : 



Q s i=Q^+Y f J, 



9e 



w b x > 



(12a) 
(12b) 
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where the hypercharge Yf a should be read from Table 1 in the same manner with 
Q£ (see, Eq. (^)). As a notable example, one can see from Table 1 that the effective 
charge of the leptons (L and e c ) disappears in the r\ model if 6 is taken to be 
1/3 0. 

Now, due to the Z-Z' mixing, the observed Z boson mass mz x at LEP1 or SLC 
is shifted from the SM Z boson mass mz'- 



Am 2 = to 



< 0. 



(13) 



The presence of the mass shift affects the T-parameter 111 at tree level. Following 
the notation of Ref. the T-parameter is expressed in terms of the effective form 
factors g%(0), gly (0) and the fine structure constant a: 



aT=l 



a (T SM + T new ) , 



(14a) 
(14b) 



where T gM and the new physics contribution T new are given by: 



olT\ 



SM 



9 2 w(0) 



• w 



_ Am 2 



(15a) 
(15b) 



It is worth noting that the sign of T new is always positive. The effects of the Z-Z' 
mixing in the Z-pole experiments have hence been parametrized by the effective 
mixing angle £ and the positive parameter T new . 

We note here that we retain the kinetic mixing term S as a part of the effective 
Z\ coupling Q^£ in Eq. (12a). As shown in Refs. |ljjl^, ^l], the kinetic mixing term 



5 can be absorbed into a further redefinition of S and T. Such re-parametrization 
may be useful if the term Yf a 5 in Eq. (12a) is much larger than the Z' charge . 



In the Eq models studied in this paper, we find no merit in absorbing the Yf5 term 
because, the remaining Q^f term is always significant. We therefore adopt 



as the effective Z\ couplings and T new accounts only for the mass shift (|13|). All 
physical consequences such as the bounds on £ and to.z 2 are of course independent 
of our choice of the parametrization. 

The two parameters T now and £ are complicated functions of the parameters of 
the effective Lagrangian (Q). In the small mixing limit, we find the following useful 
expressions 



cyT nr;w — 



9e m z 
9z m z > 
gE rn z 

gz mz> 



1 + 0&) 



1 + 



(16a) 
(16b) 
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where we introduced an effective mixing parameter £ 

< = 9_zrn^ L _ s (17) 
9E m% 

The Z-Z' mixing effect disappears at £ = 0. Stringent limits on mz> and hence 
on mz 2 can be obtained through the mixing effect if £ is O(l). We will show in 
Sec. 5 that £ is calculable once the particle spectrum of the model is specified. The 
parameter £ plays an essential role in the analysis of Z' models. 

In the low-energy neutral current processes, effects of the exchange of the heavier 
mass eigenstate Z<i can be detected. In the small £ limit, they constrain the contact 
term 5|/c 2 m| 2 . 

3. Electroweak observables in the Z' model 

In this section, we briefly discuss the theoretical framework ill to calculate the 
electroweak observables which are used in our analysis. The experimental data 
of the Z-po\e experiments, the W-boson mass measurement and the low-energy 
experiments used in this paper are summarized in Table |^. 



The pseudo-obser 
effective coupling g£ 



Dies of the Z-pole experiments are expressed in terms of the 
, where / denotes all the SM fermions except for the top- 
quark, and a being their chirality, L or R. Following our parametrization of the 
Z-Z 1 mixing (10a), the effective coupling g* a in the Z' models can be expressed as 



g f a = (glhM + Q f £t (18) 

The SM prediction for the effective coupling (g[,)sM can be expanded in terms of 
the gauge boson propagator corrections Ag z and As 2 : 

(g f a )sM = a + bAg 2 +cAs 2 , (19) 

where the numerical coefficients a, b and c are given in Refs. 0,0. Two parameters 
Ag z and As 2 in Eq. ( [l9| ) are defined as the shift in the effective couplings g 2 z {m 2 z ) 
and s 2 {m 2 z ) E3 from their SM reference values at m t — 175 GeV jjid m H = 
100 GeV. They can be expressed in terms of the S and T parameters t3 as 

Ag% = f z (m 2 Zl ) - 0.55635 = 0.00412AT + 0.00005[1 - (100 GeV/m H ) 2 ], (20a) 
As 2 = s 2 (m 2 Zi ) - 0.23035 = O-OOSeOAS - 0.00241AT - 0.00023x Q , (20b) 

where the expansion parameter x a is introduced to estimate the uncertainty of the 
hadronic contribution to the QED coupling l/a(m 2 Zi ) = 128.75 ± 0.09 0: 

l/a(m 2 ) - 128.75 

x„= Zl . (21) 

0.09 v ; 

Here, AS, AT, AU parameters are also measured from their SM reference values 
and they are given as the sum of the SM and the new physics contributions 

AS = AS SM + S ncw , AT = AT SM + T ncw , AU = AU SM + U acw . (22) 
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Table 2. Electroweak measurements at LEP, SLC, Tevatron and LENC experiments. Tl*e 
average VK-boson mass is found in Ref. Except for the weak charge of cesium atom u, 
the data of nxpoi jpiipits given in this table, which have been reduced from the 

original data EJ'E3E3'E3Sc3'l3, is summarized in Refs. The pull factors are given at 
the best fit point of the SM and four Z' models at m H — 100 GeV and a constraint 
Tncw > 0. The best fit values of parameters in both the SM and Z' models are shown 
in Table ^| Correlation matrix elements of the Z line-shape parameters and those for 
the heavy-quark parameters are found in Ref. The data R e and are obtained by 
assuming the e-fi-r universality. 





data 


SM 


X 




V 




Z parameters @ 














m 


z (GeV) 


91.1867 ± 0.0021 












r 


z (GeV) 


2.4939 ±0.0024 


-1.0 


-1.1 


-1.0 


-0.9 


-1.0 




^(nb) 


41.491 ±0.058 


0.4 


0.4 


0.6 


0.4 


0.5 




Ri 


20.765 ±0.026 


0.5 


0.5 


0.3 


0.5 


0.5 




A 0,£ 
^FB 


0.01683 ± 0.00096 


0.6 


0.6 


0.6 


0.6 


0.6 




A T 


0.1431 ± 0.0045 


-0.9 


-0.9 


-0.9 


-0.9 


-0.9 




A e 


0.1479 ± 0.0051 


0.1 


0.2 


0.1 


0.1 


0.2 




R b 


0.21656 ± 0.00074 


1.0 


1.0 


1.1 


1.1 


1.1 




R c 


0.1735 ± 0.0044 


0.3 


0.3 


0.3 


0.3 


0.3 




A o,b 

FB 


0.0990 ± 0.0021 


-2.0 


-1.9 


-2.1 


-2.0 


-2.0 




A 0,c 
FB 


0.0709 ± 0.0044 


-0.7 


-0.6 


-0.7 


-0.7 


-0.6 




A°LR 


0.1510 ± 0.0025 


1.5 


1.6 


1.5 


1.5 


1.6 




A b 


0.867 ±0.035 


-1.9 


-1.9 


-1.9 


-1.9 


-1.9 




A c 


0.647 ±0.040 


-0.5 


-0.5 


-0.5 


-0.5 


-0.5 


W mass El 
















m i 




80.410 ± 0.044 


0.8 


0.8 


0.8 


0.8 


0.8 


LENC exp. t 


















^SLAC 


0.80 ± 0.058 


1.0 


1.0 


1.0 


1.1 


1.0 




^CERN 


-1.57 ± 0.38 


-0.4 


-0.4 


-0.4 


-0.4 


-0.4 




^Batcs 


-0.137 ± 0.033 


0.5 


0.4 


0.5 


0.4 


0.4 






-0.94 ± 0.19 


-0.3 


-0.3 


-0.3 


-0.4 


-0.3 


Qw{ 




-72.06 ± 0.44 


2.2 


-0.1 


2.2 


0.1 


0.0 






0.3247 ± 0.0040 


-1.5 


-1.4 


-1.5 


-1.5 


-1.4 




KcCFR 


0.5820 ± 0.0049 


-0.5 


-0.4 


-0.4 


-0.5 


-0.4 




9ll 
9lr 


-0.269 ± 0.011 


0.4 


0.1 


0.1 


0.4 


0.1 




0.234 ± 0.011 


0.1 


0.0 


0.4 


0.2 


0.2 


A. mi n 




23.8 


18.3 


23.5 


19.2 


18.4 



Looking for Z' bosons in Supper symmetric E§ Models through Electroweak Precision Data 9 



A convenient parametrization of ASgM, ATsm and A£/g M in terms of m t and m H 
has been given in Ref. |32|. The formulae of the Z-pole observables listed in Table || 
in terms of g^ can be found in Refs. [|^|. 

The theoretical prediction of m w is given as E3€3 

m w { GeV) = 80.402 — 0.288 AS + 0.418 AT + 0.337 AU + 0.012 x a , (23) 

by using the same parameters, AS, AT, AU (p2|) and x a (^i"|). 

The observables in the LENC experiments which are used in our analysis are as 
follows - (i) polarization asymmetry of the charged lepton scattering off nucleus tar- 
get, (ii) parity violation in cesium atom, (iii) inelastic ^-scattering off nucleus target 
and (iv) neutrino-electron scattering. Theoretical expressions for the observables of 
(i) and (iil are conveniently given in terms of the model-independent parameters 
C\ q , C2 q E2l and C% q H. The ^-scattering data (iii) and (iv) are expressed in terms 
of the parameters g^J ■ In the Z' models, these model- independent parameters can 
be written as follows: 

C tq = (C tq ) SM + AC iqi (24a) 
& f = (9Z f )sM+AgZ f , (24b) 

where the first term in each equation is the SM contribution which is parametrized 
conveniently by AS and AT The second terms ACi q and Ag^J represent the 
additional contributions from the Z-Z' mixing and the direct Zi exchange, which 
are proportional to £ and g E J 'c^m| 2 , respectively. The theoretical prediction of the 
LENC observables in terms of ACi q and Ag^J can be found in Ref. 0. 

4. Constraints on Z' bosons from electroweak experiments 



4.1. -^-analysis on the Z' models 

There are six free parameters in the Z' models - the tree level contribution to the 
T parameter T now , the Z-Z' mass mixing angle £, the direct Z2-boson contribution 
to the low-energy processes g E / 'c^,m| 2 , and the three SM parameters,./?^, a s {jnz^) 
and a{m Zi ). Throughout our analysis, we use m t = 1 73.8±5.2 (GeV) E9, a s (m Zi ) = 
0.119 ± 0.002 0, and l/a(m%J = 128.75 ± 0.09 as constraints on the Sm' pa- 
rameters. The Higgs mass dependence of the results are parametrized by x H = 
ln(m H /100 GeV) in the range 90 GeV < m H < 150 GeV. The lower bound is 
obtained at the LEP2 experiment Ej. The upper bound is the theoretical limit on 
the lightest Higgs boson mass in any supersymmetric models that accommodate 
perturbative unification of the gauge couplings 11. 

We summarize the results of the fit for the ip, Xi V an d v models: 

(1) X model (6 = 0) 
T ncw = -0.059 + 0.Ux H ± 0.098 

£(10~ 4 ) = 0.02 - 9x H ± 4.05 \ ,xU a = 17.8 + 1.2a; H , (25a) 

9%/dm% = 0.237 + 0.0032x ff ±0.107 
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Table 3. Summary of the best fit values of the parameters in the SM and four Z' models 
(S — 0) for m H = 100 GeV and Tnew > 0. The mixing parameter £ and the contact term 
g%/c^ c 'm% 2 are gi ven i n units of 10 -4 and TeV -2 , respectively. The constraints on m t and 
a a (mzi) are found in Ref. ^ while that on ct(m 2 Zl ) is given in Ref. 



Parameters 


Constraints 


SM 


X 


1> 


V 


V 


m t (GeV) 


173.8 ± 5.2 


170.5 


171.3 


170.8 


170.9 


171.3 


a s {m Zl ) 


0.119 ± 0.002 


0.1186 


0.1185 


0.1184 


0.1186 


0.1185 


l/o(m|J 


128.75 ± 0.09 


128.76 


128.74 


128.75 


128.73 


128.74 


T 

± new 


















£(1(T 4 ) 






0.32 


1.63 


-2.66 


0.68 








0.245 


1.41 


-0.839 


0.635 



(2) V model (6 = 0) 



T 

-L new 

£(10- 
9 2 e/c- 



-0.075 
1.2 - 1.3a 
1.31 - 



0.Ux H ±0.097 



H 



±4.; 



0.16x ff ±2.97 



= 22.9 + 1.4a; 



Hi 



(25b) 



(3) rj model (6 = 0) 



£"(io- 4 ) 

9%/clm. 



-0.062 
-2.7-6.7/ 
-0.814 



0.14x^ + 0.097 



H 



±9.4 



0.089x H ± 0.449 



1 Xmin 



18.7 + 0.7x 



Hi 



(25c) 



(4) v model (6 = 0) 

T now = -0.057 

C(10" 4 ) =-0.42 + 
9 2 E /clml 2 =0.619 



0.14x^ + 0.098 



-0.6x^ + 3.9 
0.024x H ± 0.275 



' Xmin 



18.0 + Tlx 



H ■ 



(25d) 



where d.o.f. = 20. The mixing angle £ and the contact term g%/ c x m Z2 are gi ven 
in units of 10 -4 and TeV -2 , respectively. The best fit value of T new falls into the 
unphysical region (T new < 0) for all Z' models even if the Higgs boson mass is its 
upper limit (~ 150 GeV). It should be noticed that T ncw and £ are consistent with 
zero in all models while gfj/c^m 2 ^ shows the deviation from zero in the l-cr level 
for Xi Wi v models. 

The best fit values of the six-parameters for m H = 100 GeV under the condition 
Tnew _ are shown in Table |3[ together with the SM best fit result at m H = 
100 GeV. Only the best fit value of giVc^m 2 ^ in the rj model is found in the 
unphysical region (gfj/c^m 2 ^ < 0). The pull factors of the electroweak observables 
at the best fit point are also shown in Table ||. We learn from the table that 
almost no improvement of the fit over the SM is found for the Z-po\e and m w 
measurements. However, the Xi Vi v models show the excellent fit to the weak charge 
of cesium atom QwQj^Cs): the pull factor is reduced from 2.2 (SM) to less than 
0.1. This may imply that more than 2-er deviation of the APV data from the SM 
prediction could b_e explained by the direct exchange of the Z2 boson in the low- 



energy processes EU. On the other hand, the ip model does not show the reduction 
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Figure 1: Contour plot of A X 2 = XminOfe, $) ~ Xmin(SM) for m H = 100 GeV. The 
mixing angle (3e for the Xj^jV an d v models are shown by vertical dotted lines. 
The step of each contour is 1. 

of the pull factor in Qwil^Cs). Since all the SM matter fields in the ip model have 
the same U(l)' charge (see Table [j]), the couplings of contact interactions are parity 
conserving, which makes the contact term useless in the fit to the APV. 
We introduce a parameter 

A* 2 ^ xLnt/^) " Xmin(SM), (26) 

to measure the goodness of the fit in the Z' models compared to the SM. We can 
see from Table || that the X:V ancl v models lead to A% 2 = — 5.5(x), —4. 6(77) and 
—hA(v\ respectively while the ip model shows no improvement of the fit, A X 2 = 
-0.3. 

In order to see the impact of kinetic mixing on the fit, we show the contour 
plot of Ax 2 from the electroweak data under the conditions r ncw , g% / c^m 2 Z2 > 
on the {Pe-,5) plane in Fig. |l|. We can see from the figure that the fit of the -q 
model at S = is rather worsen (A% 2 ~ —3) as compared to that given in Table | 
(A% 2 = —4.6) because the fit in the table has been found without imposing the 
constraint gfj/c^m?^ > 0. The leptophobic r\ model (6 = 1/3) does not improve 
the fit because, due to the leptophobity, the model does not have the contact term 
g%/c^rn? Z2 which is used to make the fit to the LENC (essentially the APV) data 
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8 




Pe 



Figure 2: Contour plot of the 95% CL lower mass limit of the Z2 boson obtained 
from the LENC experiments for gs = gy and m H = 100 GeV. The vertical dotted 
lines correspond to the Xi "^i V and v models. The limits are given in unit of GeV. 

better. We find that the Z' model with ((3e,8) « (— 7r/4,0.2) shows the most 
excellent fit over the SM where A% 2 < — 7. 

4.2. Lower mass bound on Z' bosons 

As we expected from the formulae for T nGW and £ in the small mixing limit 
(|l6|), the Z2 mass is unbounded from the Z-pole data at £ = 0. For models with 
very small £, the lower bound of the heavier mass eigenstate Z2 in the Z' models, 
therefore, comes from the LENC experiments. In Fig. |[ we show the contour plot 
of the 95% CL lower mass limit of Z2 boson from the LENC experiments on the 
(Pe, S) plane by setting gs = gy and m H = 100 GeV under the condition mz 2 > 0. 
In practice, we obtain the 95% CL lower limit of the Z2 boson mass TO95 in the 
following way: 



0.05 = 




where we assume that the probability density function P(g E /m 2 z ) is proportional 
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Table 4. The 95% CL lower bound of mz 2 (GeV) in the \, tp, n and i^nodels (<5 = 0) for gE = gy 
and m H = 100 GeV. Those in Ref. pi and the result of direct search E_| are shown for comparison. 





X 


il> 


7] 


V 


our results 


554 


137 


619 


342 


results in RefJ3| 


451 


136 


317 


284 


direct search 


595 


590 


620 





to exp(-x 2 (.9i/«4 2 )/2). 

We can read off from Fig. |^ that the lower mass bound of the Zi boson in 
the ip model at S — is much weaker than those of the other Z' models. This is 
because, as we mentioned before, the U(l)' charge assignment on the SM matter 
fields in the model makes the constraint from the APV measurement useless. We 
also find in Fig. || that the lower mass bound of the Zi ipson disappears near the 
leptophobic remodel (@e = tan _1 (y / 5~73) and S = 1/3) 113. Furthermore the lower 
mass bound tend to be small at the "best fit" point which we found in Fig. |l|, 
09 E ,5) = (-tt/4, 0.2). 

We summarize the 95% CL lower bound on mz 2 for the x, tp, r\ and v models 
(5 = 0) in Table ||. For comparison, those in Ref. || are given in the same table. 
It should be noticed that the bounds on Z x , Z n and Z v masses are more severely 
constrained as compared to Ref. || due to the improved value of Qwil^Cs) while 
the bound on the Z^ mass is almost unchanged. 

We have found that the .Z-pole, m w and the LENC data constrain (T now ,£), 
T ncw and g 2 E / c 2 /m 2 Z2 , respectively. We can see from Eq. ( [i"6| ) that, for a given £, 
constraints on T new ,£ and ffs/c^m^ can be interpreted as the bound on mz 2 - We 
show the 95% CL lower mass bound of the Z2 boson for m H = 100 GeV in four Z' 
models as a function of The bound is again found under the condition mz 2 > 0. 
Results are shown in Fig. ||(a) ~ ||(d) for the x,^),r\,v models, respectively. The 
lower bound from the Z-pole and m w data, and that from the LENC data are 
separately plotted in the same figure. Shown in the figure is the lower bound of 
m z 2 9Y I 9e for 9e — 9y- The bound on m^gy /gE is approximately independent 
of gE for 9eI Qy = 0.5 ~ 2.0 in each model 13. The Z2 mass is unbounded from the 
Z-pole data at £ = because the data constrain T ncw and £ which are proportional 
to C 2 and C, respectively. Then, the lower bound on mz 2 at very small £ is obtained 
from the LENC experiments and the direct search experiment at Tevatron. For 
comparison, we plot the 95% CL lower bound on mz 2 obtained from the direct 
search experiment in Fig. || In the direct search experiment, the Z' decays into 
the exotic particles, e.g., the decays into the light right-handed neutrinos which 
are expected for some models, are not taken into account. We summarize the 95% 
CL lower bound on mz 2 for the X^iV an( i v models (5 = 0) obtained from the 
low-energy data and the direct search experiment in Table |i[ The lower bound 
of mz 2 in the rj model from the LENC experiments is competitive the bound from 
the direct search experiment. 

The lower bound of mz 2 is affected by the Higgs boson mass through the T 
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-2-10 1 2-2-10 1 2 

c ? 

Fig. 3. The 95% CL lower mass limit of Z2 in the x, 1/), 77 and v models for m H = 100 GeV. 
The Z2 boson mass is normalized by g^/gy. Constraints from Z-pole experiments and LENC 
experiments are separately shown. The results of the direct search at Tevatron £3 for the x> 4> an d 
n models are also shown. 



parameter. As we mentioned previously, T ncw tends to be in the physical region 
(Tncw > 0) for large m H [x H ). Then, we find that the large Higgs boson mass 
decreases the lower bound of m^ 2 . For £ = 1, the lower mz 2 bound in the x^i 11 
(77) models for m H = 150 GeV is weaker than that for m H = 100 GeV about 
7% (11%). On the other hand, the Higgs boson with m H — 80 GeV makes the 
lower mz 2 bound in all the Z' models severe about 5% as compared to the case 
for m H — 100 GeV. Because T now and £ are proportional to £ 2 and C, respectively 
(see Eq. (fL6|)), and it is unbounded at |C| — 0, the lower bound of mz 2 may be 
independent of m H in the small |£| region. The m ff -dependence of the lower mass 
bound obtained from the LENC data is safely negligible. 

It has been discussed that the presence of Z2 boson whose mass is much heavier 
than the SM Z boson mass, say 1 TeV, may lead to a find-tuning problem to stabilize 
the electroweak scale against the U(l)' scale E3. The Zi boson with mz 2 < 1 TeV 
for qe — gy is allowed by the electroweak data only if £ satisfies the following 
condition: 

—0.5 <( < + 0.4 for the x> i>-, v models, , , 

-0.6 < C < + 0.6 for the 77 model. ^ > 



5. Light Z' boson in minimal SUSY i?6-models 
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Table 5. Coefficients of the 1-loop /3-functions for the gauge couplings in the MSSM and the 
minimal Eg models. The model x(16) has three generations of 16 and a pair 2 + 2. The model 
x(27) and tp, n, v have three generations of 27 and a pair 2 + 2. 





MSSM 


X(16) 


X(27) 




V 


V 




33 


33 


48 


48 


48 


48 


h 


5 


5 


5 


5 


5 


5 




1 


1 


4 


4 


4 


4 


h 


-3 


-3 


















6 + — 


2 

9 + fo 


9+^ 


9+ fa 2 


9+ fa 2 


biE 












-h 



It is known that the gauge couplings are not unified in the E§ models with 
three generations of 27. In order to guarantee the gauge coupling unification, a 
pair of weak-doublets, iZL.and H', should be added into the particle spectrum 
at the electroweak scale E3. They could be taken from 27 + 27 or the adjoint 
representation 78. The U(l)' charges of the additional weak doublets should have 
the same magnitude and opposite sign, a and —a, to cancel the U(l) anomaly. In 
addition, a pair of the complete SU(5) multiplet^udi as 5 + 5 can be added without 
spoiling the unification of the gauge couplings E_lE3. 

The minimal Eq model which have three generations of 27 and a pair 2 + 2 
depends in principle on the three cases; H' has the same quantum number as L or 
Hd of 27, or H u of 27. In the following we represent the hypercharge and the U(l)' 
quantum numbers of the additional pair as (—1/2, a) for H' and (+1/2, —a) for H', 
where the U(l)' quantum number a in each Z' model follows the same normalization 
in Table 0. 

In the minimal model, the following eight scalar-doublets can develop VEV to 
give the mass terms m 2 z and rn 2 zz , in Eq. (||): three generations of H u ,Hd, and an 
extra pair, H' and H' . We express their VEVs as follows: 

i=l i=l 

where i is the generation index. The sum of these VEVs gives the observed /i-decay 
constant: v 2 +v 2 + v 2 H , + Vjp- = v 2 = w (246 GeV) 2 . By further introducing 

the notation 



tan/3 = x 2 = H „ H ' , (30) 



Vd 

we can express ( in Eq. ([l7]) as i 

C = 2 j-Q^ (1 - x 2 ) sin 2 (5 + Q% d (l - x 2 ) cos 2 (3 + Qf '.t 2 | - <J. (31) 

Because if' and if' are taken from 27 + 27, the U(l)' charge of H' , , is identified 
with that of L, Hd or H„. 
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Table 6. Predictions for gg and 8 at fi = mz 1 in the minimal models and the »?bkm model \_ 
The \J(1) Y gauge coupling gy is fixed as gy = 0.36. 



model 


a 


9E 


9e/9y 


8 


X(16) 


3 


0.353 


0.989 


0.066 






-2 


0.361 


1.010 


-0.044 


X(27) 


3 


0.353 


0.989 


0.066 






-2 


0.361 


1.010 


-0.044 


V> 




1 


0.364 


1.020 


0.028 






2 


0.356 


0.999 


0.056 






-2 


0.356 


0.999 


-0.056 


V 




1/6 


0.366 


1.025 


0.018 






-2/3 


0.351 


0.982 


-0.071 


V 






0.361 


1.010 


0.044 








0.353 


0.989 


-0.066 


l?BKM 


12 




0.308 


0.862 


0.286 



Let us remind the reader that, in the x model, three generations of the matter 
fields 16 and a pair of Higgs doublets make the model anomaly free. In this case, 
£ is found to be independent of tan p: 

C = 2Q^ - 6. (32) 

We can now examine the kinetic mixing parameter S in each model. The boundary 
condition of 5 at the GUT scale is 6 = 0. The non-zero kinetic mixing term can 
arise at low-energy scale through the following RGEs: 

haj, (33a) 

(b E + 2b 1E S + b x 6 2 )al (33b) 

{biE + h5)a>i, (33c) 

where i = 1,2,3 and t = In /i. We define a\ and as 

ol x = - — , a A = -^-. 34) 

3 4tt' 3 4tt y ' 

The coefficients of the /3-functions for a%, and S are: 

6 1 = |Tr(r 2 ), 6 E = |Tr(Q|), b 1E = ^Tr(YQ E ). (35) 



d 


1 


dt a ^ 


2^ 


d 


1 


dt ai = 


2^ 




1 






dt 


2^ 



From Eq. (33c), we can clearly see that the non-zero 6 is generated at the weak 
scale if b\ E ^ holds. In Tabic |^, wc list bi,b E and b 1E in the minimal XjV'j 7 ? 
and v models. As explained above, the x(16) model has three generations of 16, 
and the x(27) model has three generations of 27. We can see from Table || that 
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Table 7: Predictions for the effective Z-Z 1 mixing parameter £ in the minimal ip, r\ 
and v models for x 2 = and 0.5, and tan/3 = 2 and 30. 







x 2 


= 


x 2 = 


= 0.5 






tan/? 


tan/3 




a 


2 


30 


2 


30 


X 


3 


-0 


.88 


0.14 




-2 




-0.77 




i> 


1 


0.60 


1.02 


0.55 


0.76 




2 


0.58 


1.00 


0.79 


1.00 




-2 


0.69 


1.11 


-0.16 


0.06 


V 


1/6 


-1.02 


-1.35 


-0.35 


-0.52 




-2/3 


-0.93 


-1.26 


-1.11 


-1.26 


V 




0.36 


0.77 


0.57 


0.77 




-V378 


0.47 


0.88 


-0.34 


-0.14 



the magnitudes of the differences b\ — &2 and &2 — &3 are common among all the 
models including the minimal supersymmetric SM (MSSM). This guarantees the 
gauge coupling unification at fi = incur — 10 16 GeV. It is straightforward to 
obtain gE{mz 1 ) and 8{mz x ) for each model. The analytical solutions of Eqs. (33a)~ 
( |33c| ) are given in Ref. ||. In our calculation, 0:3 (m^) = 0.118 and a(mx 1 ) = 
e 2 (to,^ ) /47r = 1/128 are used as example. These numbers give gy(mz x ) — 0.357. 
We summarize the predictions for qe and S at fi = mz 1 in the minimal E e models in 
Table ^. In all the minimal models, the ratio gE/gy is approximately unity and \b\ 
is smaller than about 0.07. Some further extra fields, therefore, may be needed to 
give 5 = 0.2 which leads to the "minimal A% 2 " when (3e — — t/4, which we found in 
Fig. [|. We also^how the result of the quasi leptophobic r\ model (?7bkm) proposed 
by Babu et al. 113 in the same table. The t/bkm has, besides three generation of 27, 
two pairs of 2 + 2 from 78 and a pair of 3 + 3 from 27 + 27 in order to achieve the 
leptophobity (5 ~ 1/3) at the weak scale through the quantum corrections. We find 
that the ?7bkm model predicts gs/ffy ~ 0.86 and S ~ 0.29, which is rather close to 
the leptophobity, S = 1/3. 

Next we estimate the parameter £ for several sets of tan/3 and x. In Table |?], 
we show the predictions for £ in the minimal Xi V an< i v models. The results 
are shown for tan/3 = 2 and 30, and x 2 = and 0.5. We find from the table that 
the parameter £ is in the range \C\ < 1 .35. It is shown in Fig. || that mz 2 gy /gE 
is approximately independent of gE/gy- Actually, we find in Table | and Table fj] 
that the predicted \8\ is smaller than about 0.1 and gE/gy is quite close to unity 
in all the minimal models. We can, therefore, read off from Fig. ^ the lower bound 
of mz 2 in the minimal models at gE = gy- In Table |^, we summarize the 95% CL 
lower mz 2 bound for the minimal Xi^Pj 7 ! an d v models which correspond to the 
predicted C in Table 0. Most of the lower mass bounds in Table || exceed 1 TeV. 
The Z^ boson with mz 2 ~ 0(1 TeV) should be explored at the future collider such 
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Table 8: Summary of the 95% CL lower bound of m^ 2 (GeV) which corresponds to 
the predicted £ in Table 
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30 
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30 
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1490 


620 




-2 




1380 




i> 
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1270 


1780 


1200 


1470 




2 


1250 


1760 


1510 


1760 




-2 


1380 


1890 


540 


370 


V 


+ 1/6 


1440 


1840 


660 


860 




-2/3 


1330 


1730 


1550 


1730 


V 


+ V /l /6 


1030 


1600 


1340 


1600 




-V378 


1200 


1730 


880 


580 



as LHC. The discovery limit of the Z' boson in the Eq models at LHC is expected 
as (in unit of GeV) 



X 


</> 


rj v 


3040 


2910 


2980 * * * 



(36) 



All the lower bounds of mz 2 listed in Table |8j are smaller than 2 TeV and they are, 
therefore, in the detectable range of LHC. But, it should be noticed that most of 
them (1 TeV <mz 2 ) may require the fine-tuning to stabilize the electroweak scale 
against the U(l)' scale E3. 

6. Summary 

In this review article, we have studied constraints on Z' bosons in the SUSY 
Eq models. Four Z' models — the %, tp, r\ and v models are studied in detail. The 
presence of the Z' boson affects the electroweak processes through the effective Z-Z' 
mass mixing angle £, a tree level contribution T new and the contact term g E / 'czm^ , 
where the latter two parameters are positive definite quantities. The Z-pole, m w 
and LENC data constrain (T new ,£), T nGW and g^/c^m^, respectively. From the 
updated electroweak data, we find that three Z' models (%, rj, v) improve the fit 
over the SM where the total \ 2 decrease about five units, owing to the excellent fit 
mainly to the improved data of parity violation in cesium atom which is expressed 
by the weak charge QwQ^Cs). The more than 2-er deviation of QwQ^Cs) from 
the SM prediction could be explained in these three Z' models. Due to its parity 
conserving property of the U(l)' charge assignment on the SM matter fields, the ip 
model does not improve the fit to the QwQ^Cs) data. The impact of the kinetic 
mixing (6 ^ 0) on the fit is also examined on the (Pe, 8) plane. The Z' model with 
(Pe,8) — {— 7r/4, 0.2) shows the most excellent fit to the data among the SUSY 
Eq models where the total \ 2 decreases by about seven units as compared to the 
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SM best fit. The 95% CL lower mass bound of the heavier mass eigenstate Z2 is 
shown as a function of the effective Z-Z' mixing parameter £ together with the 
result of direct search experiment. By assuming the minimal particle content of 
the E 6 model, we have found the theoretical predictions for (. It is shown that 
the Eq models with minimal particle content which is consistent with the gauge 
coupling unification predict the non-zero kinetic mixing term 5 and the effective 
mixing parameter £ of order one. The present electroweak experiments lead to the 
lower mass bound of order 1 TeV or larger for those models. 
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